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Abstract. If a single particle obeys non-relativistic QM in R** and has the 
Hamiltonian H = — A -f /(r), where /(r) = ^*L]^ n^r'', 2 < qi < qt+i, at > 0, 
then the eigenvalues E = (X) are given approximately by the semi-classical 
expression E = min s -I- y^*L, o-i (Pir)''^ > . It is proved that this formula yields 

r>Q I ^ ^1-1 J 

a lower bound if Pi = P^'g'iqi), an upper bound if Pi = P^'^' (qk) and a general 

approximation formula if Pi = P^^ (gi )- For the quantum anharmonic oscillator 

/(r) = + \r^"^,m = 2,3,... in d dimension, for example, E = i?^'^'(A) is 

determined by the algebraic expression A = i (^^j^rx^)™ ( (^e"-^) ~ (m-i) ) 

where 5 = E'^m? — ioi{m? — 1) and a,/3 are constants. An improved lower 
bound to the lowest eigenvalue in each angular-momentum subspace is also 
provided. A comparison with the recent results of Bhattacharya et al (Phys. 
Lett. A, 244 (1998) 9) and Dasgupta et al (J. Phys. A: Math. Theor., 40 (2007) 
773) is discussed. 
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1. Introduction and main results 



The purpose of the present work is to estabhsh a global bound formula for the discrete 
spectrum {E^f^^}, n = 1,2, ... J — 0,1,2, ... oi the d-dimension Schrodinger equation 
with polynomial potentials given by 

/ k \ 



-A + ^ a^r«0 V = i^^A, '2.<q^<q^+l, 



(1) 



where A is the d-dimensional Laplacian operator, r = ||r|| , r S R'', and the coefficients 
o-i > 0, are not all zero. The key motivation for our present study lies in the well-known 
fact that the majority of quantitative predictions of Schrodinger's equation with a 
polynomial potential ([1]) in nuclear, atomic, molecular, and condensed matter physics 
must usually rely on numerical estimates [I]-[S]. Thus, a simple global eigenvalue 
formula can serve as a basis for exploration and also for checking different approximate 
methods in quantum mechanics [B]. Another important motivation for the present 
work is a recent contribution by Dasgupta et al regarding a general simple scheme 
for evaluating the ground state as well the excited-state energies for Ar^™ quantum 
anharmonic oscillators in one dimension, see also [B] . We provide in the present work 
a more general scheme sufficient to generate all energy levels in arbitrary dimension, 
not only of r^™ quantum anharmonic oscillators, but also for any polynomial potential 
of the form X]i=i '^i''"'^' ^ sufficient degree of accuracy to be interesting. The purpose 
is not merely to obtain accurate energy eigenvalues for different polynomial potentials 
for which a large number of methods exist in the literature. Rather, we propose 
a simple approach which provides energy bounds as well as an approximate energy 
formula with a reasonable accuracy and with a minimum amount of effort. Consider, 
as an example, the celebrated quantum anharmonic oscillator [5]-[23j Hamiltonian 



-A- 



' Ar 



,2m 



2, 3, ... in d dimensions: we show that for any state n = 1, 2, . 



the eigenenergy E ~ E^^ (A) is determined approximately by the expression 



2a{m - 1) 
mE — 5 



4a 



E 



{mE — S) (m — 1) 



(2) 



where S — ^ E'^rrv^ — 4:a{m'^ — 1) and a and /3 are constants. Further, we show that 
upper or lower bounds for the energy eigenvalues ([2]) for a given state are expressed in 
terms of a single constant for any value of A. The dependence of a and (3 on m and d 
will be discussed in a subsequent section. We obtain our global eigenvalue formula for 
((D) by using the so-called P-representation [24] for the Schrodinger spectra generated 
by the pure power-law potential {q > 0). In this representation, a discrete eigenvalue 
e is written as the minimum of a function of one variable r and a parameter P: this 
induces a one-one relation between e and P. More specifically, we write 

(-A + r'')V^i^; = ^ ei'^q) = min ' 



r>0 



where 



(2+q)/2q 



-,1/2 



2 + q 



,(3) 



(4) 



This may seem at first sight rather inconvenient since the computation of P requires 
knowledge of e. An important advantage of (jl]), however, is that the computation of 
P is independent of the potential parameters. In other words, the computation of 
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P for _ff = — A + r"^ is sufficient to yield the discrete spectrum of the Hamihonian 
= — A + vr'' with eigenvalues given by 

, 9 



^ife il) = mill 

"'^ r>0 



^-0 



(5) 



for arbitrary v > 0. This may seem unnecessary for a Hamiltonian of the form 
because a simple scaling argument shows E^f^^ {q; v) — w^/ E^^ (g) ; but for 
polynomial potentials, as in ([T|), where ai > 0, i = 1, 2, . . . , fc, equations ©-(HI) play 
an important role in establishing some of the general energy formulae [25j through the 
decomposition of the Hamiltonian ([T]) by means of 

k k 

H = -I^ + ^a,ri^ ^^uj,H^^ (6) 



i=l 



where 



(7) 



and is an arbitrary set of positive weights with sum equal to 1. Further, it 

worth mentioning that this dependence can be resolved for certain special values of 9, 
for example if g = 2, we know that [26j 

P^^(2) = (2n + Z + f -2) ifd>2 

(8) 



0(2) 



P„(2) = (n-i) 



if d = 1 



The main results of the present work may be summarized by the following two 
theorems: 

Theorem A: Eigenvalue bounds for the spectrum {E^^^} of the Hamiltonian (QJ) are 
given by 



£ = min 

r>0 



+ ^a,(P,r)«' 
i=i 

P^fiq,), then£<E^^ 



(9) 



P(?fe) then£>E^^^ 



where 
i) ifP, 

a) if p. 

Here the numbers pI^i {qi) are given by 

Theorem B: The eigenvalues of the Hamiltonian (Qp are given approximately by the 
semiclassical formula 

r k 

- + Y.a,{P,rr ■ (10) 



£ 



mm 

r>0 



where, for the lowest eigenvalue in each angular-momentum subspace, and d > 1, we 
have 



i) For n = I and Pi = P^f\qi),', 



1 . . . fc, then £ < E^) for all d > 1 and 



I = 0,1,2, 
Forn>2 

Further, for the lowest eigenvalue in d> 1, we have 



ii) Forn>2 and P, ^ P^f{q,),i ^ 1 . . .k, £ ^ E^f . 
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Hi) £ < if the numbers Pi are replaced by the explicit lower approximations for 



Pwi^i) g^ven by 



P, 



d 



r[i + f) 



exp(l) 



(11) 



iv) £ > e''^ if the numbers Pi are replaced by the explicit upper approximations to 
Pioili) d^ven by 



(12) 



The difference between the two parts of Theorem B is that, in the first part (i)-(ii), the 
P-numbers are to be computed from the pure-power energies by use of (j4]), whereas, in 
the second part (iii)-(iv), the P-numbers are given exphcitly in terms of the Gamma 
function. We use the term "semiclassical" in the following sense: once the component 
kinetic potentials have been fixed by the P-numbers, what remains is a minimization 
over a real function; in the approximation, this expresses the trade off between the 
kinetic and potential energies; the final picture is semiclassical since the kinetic energy 
is reduced to 1 /r^ and a wave equation is no longer involved. 



In the next section, we discuss the proof of these two theorems. In section (3), the 
application of these two theorems to the quantum anharmonic oscillator Hamiltonian 
is presented. The conclusion is given in section 4. 



2. Proof of theorems A and B 



2.1. Proof of Theorem A 

The proof of theorem A depends on the application of envelope theory and Kinetic 
potentials technique developed earlier by Hall [22]- [29] and used successfully since 
then. We shall outline here a brief summary of the theory to provide us with sufficient 
details to prove the theorem, and we refer the interested reader to Ref. [57]- [29] for 
more details. For simplicity, we present this brief summary for the case of d = 3 
spatial dimensions: for arbitrary d, the extension is straightforward. Consider the 
Schrodinger operators of the form 

H = -A + vf{r), (13) 

where / is the shape of a central potential in R'^, and u > is the coupling parameter. 
The principal idea of envelope theory [22]- [2H] is that the minimization of the Rayleigh 
quotient {ip, H^p) / {ip, tp) is performed in two stages. The first stage, with {tp, —Atp) = s 
fixed, involves only the shape of the potential / and leads to a family {/„£} of kinetic 
potentials /„^(s). Here s is a positive constraint variable: it only becomes the mean 
kinetic energy when the minimization of the sum of the kinetic and potential energies 
has been effected. We have 

Ene = min{s + vj^(,{s)} (14) 

s>0 

in which the critical value of s — {ip, —Aip) > is the mean kinetic energy. The kinetic 
potentials [25j . which represent the result of min-max theory applied to the potential 
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shape / for fixed s, are given as a Legendre transformation {s — E{v) — vE'{v), f{s) — 
E'{v)} of the function E{v), which describes how the eigenvalue depends on the 
coupling V. They may also be defined by the following general formula 

7„,(s) = inf sup /,AM/([(^,-A^)/s]i/2^)V,(r)dV, (15) 

\\ti>\\=i 

where Vni is the span of a set of n linearly independent functions. It is interesting 
to notice that the kinetic potential /„£ can be replaced by the potential /(r) itself 
through the parameterization of f ni{s) in term of the variable r (used here as a new 
parameter to replace s), that is to say fnt{s) — f{r). We now invert this monotone 
function to give the K functions [IHj 

s-(7:/°/)m-<^'m- (16) 

It is easy to show that the K functions obey the scaling property 

and in general they are independent of coupling and potential shifts 1291 . The 
eigenvalues are recovered from the K functions by the expression 

E.,u - Fnt{v) = min {i^^^^ (r) + vf{r)] (18) 

For the power-law potentials /(r) = r^, it is known by mean of simple scaling argument 
that the spectrum of the pure-power Hamiltonian satisfies 

_ A -f z;r« ^ fI^^{v) = E^S{iy^'i+^\ (19) 

In order to compute the kinetic potentials fnds), we notice from the minimization 
process of (|14p that fni{s) = —v^^ and consequently we have 

. = Fifiv) - v7^,is) ^ 7„,(5) = ^F^fiv) (20) 
which implies using that 

7n,is) = ^v-^El^,\l). (21) 
On the other hand, we have from the l.h.s. of (l20l) that 



v''Fifi^)^fni{s)~sfUs) (22) 



which implies using ([2T|) that 



The K functions are then computed by means of and (|20p - (P5|) to yields 

(n-(i^V'™''fM^.)" 



where we have defined 



Pni{q) = (i?, 



(9) 



(2+g)/2g 



1/9 



1/2 



(25) 
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The eigenvalues are then recovered by (|18p as 

In ordered to obtained a definite bound, Hall [57] used interesting geometric 
interpretation in terms of envelopes. If the potential shape /(r) = g{h{r)) is a smooth 
transformation g of a soluble potential h, then the kinetic potentials associated with 
/(r) are given by 

/(r) = g{h{r)) ^ 7„,(s) « giK.is)) (27) 

and the corresponding K functions satisfies 

Kl!f^ - (5 ° h^ir' o{goh) = kI o = (28) 

Therefore 

/ g{h) ^ K'- f^ = i^C^) (29) 
and the eigenvalue approximations are given by 

E„,^mm{Kl',\r)+vf{r)} (30) 

in which g no longer appears. This expression yields upper or lower bounds depending, 
respectively, whether g is concave or convex [25]-[26]. For /(r) = g{r'') = X^iLi o.ir'^': 
since qi < Qi+i, clearly g is convex if q — gi (lower bound) and concave if g = qk 
(upper bound). We therefore have, by using (j30p with h{r) = r'^, 

^;„..min|(^)%.^a..^j. (31) 

Or, equivalently, and by a change in the minimization variable, 

E^t = Tam\^ + vy^a,{Pni{q)rY'\ . (32) 

I i=l ) 

With V = 1, this equation yields Q and we obtain a lower bound if g = gi and an 
upper bound if g = g^ . This complete the proof of the theorem. D 



2. 2. Proof of Theorem B 

The first part of theorem B was introduced [SS] to improve the lower bounds for the 
ground state energy obtained in theorem A. The second part is based on the Barnes 
et al's [30] general lower-bound formula for the lowest eigenvalue of the Schrodinger 
operator H = —A + V{r) in d > 1 spatial dimensions. The extension to the potential 
sums, such as that of ([T]), was introduced in 31j where a detailed proof of theorem B 
can be found. □ 
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3. Fractional anharmonic oscillator 



Before we study specific problems in quantum mechanics, we first consider the 
application of theorem A and Theorem B to the class of arbitrary fractionally 
anharmonic oscillator Hamiltonians |32j : 

H = -A + Y, 9sr' (33) 
Sez 

where Z is an arbitrary finite set of the integer or rational numbers and the coupling 
gs,5 e Z are chosen so that the Hamiltonian supports the existence of a discrete 
spectrum. It is known |32j-[35j that this class of Hamiltonian possesses elementary 
solutions for certain particular cases of the coupling gs. For consistency, we assume 
6 > 2, although the conclusion of theorems A and B are perfectly applicable for 
all S > —1, where, for example, the P number in the case (5 = — 1 is P!^^{—^) = 
(n + / + d/2 — 3/2). This class of Hamiltonian is a generalization of the Hamiltonian 

= -A + Cr" + Dr^, /3 > a > (34) 

which has been used in the theory of heavy quarkonia |36|-[37|. Denote q = min{(5} 

and Q — max{(5}. By using theorem A, we immediately find analytic expressions for 

lower bounds ej^^ and upper bounds Ef^^ for the eigenvalues of the Hamiltonian ((33|) : 
these can be written explicitly as 



r>0 



-+ ^ gmq)rf 



and 



mm 

r>0 



5GZ 
5=min{5} 



E 

Q=max{5} 

5ez 



(35) 



9s{P{Q)r) 



(36) 



Here the numbers P{q) and P{Q) are computed numerically by means of Eq.Q for 
rational g, Q ^ —1,2 by the use of direct numerical integration of the corresponding 
Schrodinger equations (—A + r'^)ip = Eqip and (—A + r'^)V' = Eq^p. An interesting 
improvement for the eigenvalues £^^((5) and Efg{S) can be obtained through the 
application of theorem B. The cost, however, is that the exact eigenvalues of the 



rational power-law potentials V{r) 



for each S € Z must be computed numerically. 



Less accurate bounds can be obtain directly using the explicit P numbers (jlip and 
p2p . An important class [35^ of the fractional anharmonic oscillator Hamiltonians 
(|33P that have found many applications ^32 in quantum field theory [33] is given by 



H = -A + V{r) = -A 



29+1 



92q+l 



= > 0. 



(37) 



This class of Hamiltonians has found many applications not only in quantum 
mechanics (where V{r) represents [32. an arbitrary potential in the limit q oo) 
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but also, ofr example, in the Reggeon field theorem on the lattice [35]. Theorem A 
gives immediate lower and upper bounds to the eigenvalues of (I37p as: 



and 



^niiS) ~ mill 

r>0 



E^S) « min 

r>0 



2q+l 



2q+l 



(38) 



(39) 



where P^ii'^) is given by Eq.® and PniC^Q + 1) is given by ([3]), respectively. 



4. Quantum anharmonic oscillator 

In this section, we consider the Schrodinger equation 

{~LuA + ar^ + br^"')iP ^ E{Lu,a,b)i;, m = 2,3,4,... (40) 

where ui, a and b are positive parameters and the potential in (|33p is a single- well 
potential which describes for m = 2, 3, 4, 5, . . . the quartic, sextic, octic, and decadic 
oscillators, and so on. It is easy to check that for the energy in PO)) the following 
scaling relation holds 

/ ^|_^(m-l)/2^ 

i?(^,a,6) =(ac.)V2i5(^i,i,___ 
^{r;.,a,b)=^j[^(-yr;l,l,-^,^y (41) 



Thus the original problem (|40p is essentially a single-parameter problem which we now 
write as 

(-A + r2 + Ar2")V' = £;(A)V', to = 2,3,4, ... (42) 

where E{X) = i?(l,l,A) and A — ''^(L+i)/2 ■ The Schrodinger equation with the 
quantum anharmonic oscillators (|42p are among the most widely studied models in 
quantum mechanics. In spite of their simplicity, they give rise to interesting problems, 
both computationally and conceptually f23' . A rigorous analysis of the mathematical 
properties of the anharmonic oscillator Hamiltonians was made by Simon [8^ and 
by the seminal work of Bender and Wu [9 . The aim in the discussion we present 
below is to derive simple upper- and lower-bound formulae based on Theorems A and 
B. For the anharmonic oscillator potentials 

/(r) = Ar^", m = 2,3,... (43) 

Theorem A implies that 

" 1 



£{\) ~ min 

r>0 



ar'^ + Xfir 



(44) 



where 

• £(A) < E{\) is an lower bound, if a = {Pne{2)f and (3 = (P„f (2))2™. 

• £{\) > E{X) is an upper bound, if a = (F„f (2m))2 and f3 = (F„f (2m))2™. 
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Furthermore, Theorem B imphes that if a = (P„f(2))2 and (3 — (P„f (2™))^'", then 
£{X) < E{X) for n = 1, and £{\) w E{\) for all n> 2. Let x — r^, we note that the 
minimization of ([33]) occurs at —-^+ct + mX(3x"''~^ = 0. Multiplying through by x 
and solving for Xj3x™' , we can easily show that the minimization of (j44p occurs at 



2 mE — ^jE'^m'^ — Aa{m'^ — 1) 



r = ^hr^ (45) 

2a{m - 1) ^ ^ 



and consequently we have 



1 /2a(m- 1)\ /4a E 



^ I3 \ mE-S ) \(mE-6) (m-1)/ ^^^^ 



where S — y^E^m'^ — Aa{m'^ — 1). Thus for finding the energy eigenvalues of 
anharmonic-oscillator Hamiltonians in ((42|) one has to solve Eq. (|46| for the given 
A. It is clear that at A = 0, equation (|46| implies E — 2y/a, with a — (pjj^(2))2 = 
(2n + / + I - 2)2, as given by Consequently, = 4n + 2Z + d - 4, the resuh for 
the d-dimensional harmonic oscillator [5S]. When equation ([^5]) is used to determine 
the lower or the upper bounds to the exact eigenvalues of the Ar^"* oscillator, it is 
clear that the formula is expressed in terms of a single constant for any value of A. 
This follows from the fact that, for lower or upper bound, a™ — (i and equation p6)) 
reduces to 



2™(m-l)(™-i) i-E + Jm'^E'^ ~ Aaim'^ - 1)) , , 

A = ^ ^ ^ ^ " . 47 

(to + 1) (jnE - ^to2£;2 - 4a(m2 - 1))™ 

This is a remarkable simple formula that gives a global lower and upper bound to the 
exact eigenvalues for a given A for all n = 1, 2, . . . and I ~ 0,1, 2,... in d dimensions, 
accordingly as a = (2n + ? + | — 2)2 and a = (P^^^(2to))2, respectively. In particular, 
a global formula that gives a lower bound for all n = 1, 2, . . ., to = 2, 3, . . . is 



2™(m-l)(™-i) + y^m2£:2 - 4(2n + Z + f - 2)2(m2 - 1)) 

A — ^ — . 

(™ + ^> {mE - ^m^E^ - 4(2n + ? + f - 2)2(to2 - 1))™ 

(48) 

Note in the case of d = 1, we should set either Z = — 1 or i = to obtain a lower 
bound to even or odd (exact) eigenvalues, respectively. Despite the generality of 
(|47p. we should like to make two immediate remarks concerning the application of 
theorem A: (i) Formula (|T7|) , in general, gives a loose bound; (ii) the upper bound a = 
(P^^^(2to))2, TO = 2,3, .. . which is obtained by means of Eq.Q, requires knowledge 
of the exact eigenvalues of Schrodinger equation (— + r'^™)ip = e'"^) {2m)%lj . In 
this paper we have found the values of P^^^ (2to) by the numerical integration of 
the Schrodinger equation just mentioned, and then we used Eq.Q to find the 
corresponding P-numbers. For immediate use of equations l|47p and (|48p . we report 
in Table 1 the values of P^''q''(2to) and (3 = (P^J''(2to))2'" for different values of to. 

Table 1: Values of Pi^o^(2m) and /3 = (P^J^ (2to))2"' for different to. 
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m 






2 


0.648 283 101 647 721 4 


0.176 627 696 530 967 9 


3 


0.752 213 287 729 753 3 


0.181 153 198 043 223 7 


4 


0.830 692 879 447 472 3 


0.226 737 786 349 046 1 


5 


0.892 746 975 167 740 8 


0.321 576 181 371 282 8 


6 


0.943 407 187 840 825 1 


0.497 038 660 113 318 



In order to illustrate the above discussion, we consider the case of finding the 
eigenvalues of the anharmonic oscillator Hamiltonian — ^ + r^ + O.Olr'*. Eq. ((48|) 

gives a lower bound 1.00248 and Eg. ([17)) with a — (P^'q' (2m))^ gives an upper bound 
of 1.32038. The exact eigenvalue in this case reads 1.00737. In order to improve these 
bounds, we can make use of Theorem B. For the ground state eigenvalues in d > 1, 
the first part of Theorem B can be applied to obtain a a more accurate lower bound 
formula. Further, the second part of theorem B can be used to obtain straightforward 
lower and upper bound without the necessity of the numerical computation of P- 
numbers, thanks to the explicit approximate values of p''^ i^m) given by ini) and 
(|12|1 . In either case, equation (j46|) gives a simple general formula for the energy bound 
of -E = i?(A) with reasonable accuracy 

1 (m-l)('"-i) (-i; + ^m2(£;2-rf2) + ^2) 



(m + 1) 



(mE 
Ml 



y'm^iE^ ^d^)+d^y 



(49) 



for m = 2, 3, . . . where now if /3 = (P-^g (2™))^™ as given by Table (1), Eg. pS)) gives a 
lower bound for given A. The results of this formula is illustrated in the last columns 
of Tables 2 and 3. On other hand, if /3 is a fixed number given by PT|) and (fT^ . then 
Eq. (|49|) gives lower and upper bound, respectively. Note that Theorem B, still allows 
us to conclude that Eq. (|49)) yields a reasonable approximation to the excited-state 
energies for n > 2. However, in this case (3 = {P^Q\2m))'^™ is strictly given by 
In the case of d = 1, Ea. (|^ reads (m = 2, 3, . . .) 



/3A, 



(50) 



1 (m - l)("-i) {~E + ^m'^{E'^ - 1) + 1) 
2^ (m + 1) {niE- ^P^{E'^ - 1) + 1)" 

where /3 = (P^^g^ (2to))2™ is given by pip and (|12p for a lower and an upper bound 
respectively. Equation ([50)1 can be compared with the recent formula introduced by 
Bhattacharya et al [5] for the approximate ground state energy of the Hamiltonian 
— A + + Ar^™ in one dimension, namely 

_ (•£i(m)^|(m-l)(l+2/(m+2+A)) ^ (if^^™) ) (™+l) A (51) 

where Kf^ = 1.06036209, Kf^ = 1.14480245, ivT^"*^ = 1.22582011 are the coefficients 
of the first terms in the respective strong coupling expansion computed by Weniger 
|23| . In Tables 2 and 3, we have compared our lower and upper bounds given by ()49p . 
for the quartic and sextic anharmonic oscillator, along with the exact eigenvalues 
obtained by the direct integration of the corresponding Schrodinger equation. We 
have also compared the best lower bound obtained using ([49)) . where /3- values 
were given by table 1, with the approximate eigenvalues of the ground state energy Eh 
computed by Bhattacharya et al [6] using formula (|50p . Recently, Dasgupta et al [7] 
have extend the work of Bhattacharya et al |S] to evaluate the excited state energies, 
still in the one-dimensional case. As in the case of the ground state ([SO)) . they found 
that the excited-state energies for the Ar^™ oscillator defined by the one-dimensional 
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Hamiltonian operator H = —d^/dr"^ + Ar^™ are also a polynomial equation of the 
same degree and are given by 



where is the n^^ excited state energy of the Ar^™ oscillator and are 

constants [7]. Our formulas and (|T7)) are more general and seems to yield more 
accurate results, even for large values of the coupling parameter A. 

5. Conclusion 

The application of envelope theory and kinetic-potential techniques to polynomial 
potentials has yielded fairly general and good energy bounds for arbitrary values of 
the coupling constants. As specific examples, the application of theorems A and B to 
the quantum anharmonic oscillators has produced a global energy formula sufficient to 
generate all energy levels in arbitrary dimension for r^™ anharmonic oscillators with a 
fair degree of accuracy. The main emphasis of this paper has been on energy formulas 
that are also bounds. However the energy formula of Theorem B (ii), namely 



which does indeed yield a lower bound for the bottom of spectrum (n = 1) in each 
angular-momentum subspace, is a remarkably general and accurate approximation: it 
requires the input of the pure-power P numbers, and it then predicts approximately, 
for all the eigenvalues in all dimensions, how the spectrum generated by the potential 
sum depends on the mixing parameters {a^}; it also has the attractive collocation 
property that it is exact whenever all but one of the potential coefficients are zero. 
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Table 2: Calculated values of upper and lower bounds, using ^3\ l, to the ground state 
energy of the quartic anharmonic oscillator, along with exact values, for different 
values of A. The comparison between the lower bound Ej^ given by Eg. using 
the exact values of P by means of table (1), and the approximate eigenvalues Eh of 
Bhattacharya et al '^61 using are also shown. 



A 


Exact value 


Lower bound 


Upper bound 


Eb 


El 






using Eq.(|43l) 


using Eq.(|43l) 


Rcf. [6^ 


using Eg. ((43)1 


0.001 


1.00075 


1.00062 


1.00075 


1.00079 


1.00071 


0.01 


1.00737 


1.00614 


1.00739 


1.00783 


1.00697 


0.1 


1.06529 


1.05585 


1.06620 


1.07005 


1.06275 


0.2 


1.11829 


1.10288 


1.12062 


1.12702 


1.11473 


1.0 


1.39235 


1.35510 


1.40332 


1.41155 


1.38754 


4.0 


1.90314 


1.83699 


1.92881 


1.91489 


1.89895 


10.0 


2.44917 


2.35648 


2.48862 


2.45005 


2.44575 


50.0 


4.00399 


3.841639 


4.078522 


3.99621 


4.00182 


100.0 


4.99942 


4.79395 


5.09516 


4.99161 


4.99766 


1000.0 


10.63979 


10.19449 


10.85151 


10.63521 


10.63896 


2000.0 


13.38844 


12.82706 


13.65591 


13.38474 


13.38778 



Table 3: Calculated values of upper and lower bound to the ground state energies 
of the sextic anharmonic oscillator along with exact values for different values of A. 
The comparison between the lower bound El given by Eq. \4-3^ and the approximate 
eigenvalues Eb of Bhattacharya et al fQ/ using ^44\ l o^so shown. 



A 


Exact value 


Lower bound 


Upper bound 


Eb 


El 






using Eq.(|43| 


using Eq.(|43| 


Ref. [6^ 


using Eq.(|43| 


0.001 


1.00185 


1.000932 


1.001859 


1.00143 


1.00144 


0.01 


1.01674 


1.008994 


1.017387 


1.01374 


1.01366 


0.1 


1.10908 


1.070681 


1.119935 


1.10565 


1.09920 


0.2 


1.17389 


1.119782 


1.192805 


1.17513 


1.16261 


1.0 


1.43653 


1.334560 


1.484050 


1.44870 


1.42400 


4.0 


1.83044 


1.675050 


1.916177 


1.83193 


1.82058 


10.0 


2.20572 


2.004582 


2.322916 


2.19235 


2.19734 


50.0 


3.15902 


2.850163 


3.348809 


3.13471 


3.15304 


100.0 


3.71698 


3.347427 


3.946987 


3.69348 


3.71187 


1000.0 


6.49235 


5.828630 


6.914382 


6.47694 


6.48941 


2000.0 


7.70174 


6.911387 


8.205757 


7.68861 


7.69925 
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